REVERSIBLE ARITHMETIC CODING FOR QUANTUM DATA 

COMPRESSION 
BACKGROUND OF THE INVENTION 



Field of the Invention 

The present invention generally relates to a method of compressing and 
decompressing a block of symbols (a block quantum state) emitted by a 
memoryless quantum Bernoulli source. This presents a simple-to-implement 
quantum system for projecting, with high probability, the block quantum state 
onto the typical subspace spanned by the leading eigenstates of the source's 
density matrix. 

Description of the Related Art 
Computers and information technology electronically transfer large 
amounts of data. The data could be stored on a computer disk or in the memory 
of a computer. The data could be transferred from computer to computer via 
networks, the internet or over the airways. Today's demands on data transfer 
range from extremely complex items such as satellite photographs, which require 
a large amount of memory to simpler items such as e-mail text messages, which 
demand less memory. 

Conventional technology performs the data transfer and storage by 
reducing the data to an electronic charge or bit which is stored in the computer 
memory, a CD ROM or disk. The bits are assigned data using a binary code, 
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which includes only two values, 0 and 1 . Thus, each bit may exist as either a 0 or 
a 1 and a series of bits may be required to represent each piece of data. 

The methodology of this approach to data transfer is derived from classical 
physics. According to these sciences, modern information theory assumes that an 
information bit can exist in only one of two states, say, 0 or L However, classical 
physics is known to fail spectacularly under many circumstances, for example, 
when the objects being described are very small or have very large energies. As a 
result, conventional information theory fails to properly describe how information 
may be represented and transformed in such physical systems. 

Quantum information theory, which is based on the study of quantum 
mechanics, is capable of more accurately describing and transferring a wider 
range of data than the conventional binary code method. Similar to the 
conventional methods, quantum information theory functions by assigning the 
data to a series of bits, but the bits are called quantum information bits. 

In contrast to the classical information bit, a quantum information bit can 
exist in a superposition of two orthogonal quantum states. Meaning 
rather than being limited to one of two states, 1 or 0, a quantum information bit 
can exist in an infinite number of states. A quantum information bit may be any 
combination of 0 and 1. For example, it can be 40% 1 and 60% 0. 

Quantum information can, in principle, provide significant advantages for 
certain problems. For example, quantum algorithms for calculating discrete 
logarithms, see, Shor, (P.W. Shor, Algorithms for quantum computation: Discrete 
logarithms and factoring " in Proceedings of the 35th Annual Symposium on 
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Foundations of Computer Science, Santa Fe, New Mexico, (Los Alamitos, CA), 
pp. 124-134, IEEE Computer Society Press, 1994, incorporated herein by 
reference), and searching unsorted databases, see, Grover, (L.K. Grover, "A fast 
quantum mechanical algorithm for database search," in Proceedings of the 28th 
Annual ACM Symposium on Theory of Computation, Philadelphia, Pennsylvania, 
pp. 212-219, 1996, incorporated herein by reference), have been discovered which 
are faster than their classical counterparts. Quantum bits, in contrast to classical 
bits, cannot be copied perfectly, and this is useful in such tasks as quantum 
cryptography, see, Bennett (C.H. Bennett, G. Brassard, and A.K. Ekert, "Quantum 
cryptography," Sci. Am., vol. 267, pp. 50-57, Oct. 1992, incorporated herein by 
reference). Furthermore, Fuchs (C. Fuchs, "Nonorthogonal quantum states 
maximize classical information capacity," Physical Review Letters, vol. 79, pp. 
1 162-1 165, 1997, incorporated herein by reference), has shown that, rather 
unexpectedly, there exist certain quantum communication channels for which the 
optimal classical information transmission rate is achieved only using 
non-orthogonal quantum states as the symbols. Finally, quite surprisingly, 
quantum error correction codes have been developed, see, Calderbank et al. (A.R. 
Calderbank, E.M. Rains, P. W. Shor, and NJ.A. Sloane, "Quantum error 
correction via codes over gf (4)," IEEE Trans. Inform. Theory, vol. 44, pp. 
1369-1387, July 1998, incorporated herein by reference) and references therein. 
Such codes might provide the key technology needed to prevent decoherence of 
quantum states, and, hence, a way to realize large-scale quantum computing 
devices. For excellent reviews of the field, see, for example, Bennett and Shor 
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(C.H. Bennett and P.W. Shor, "Quantum information theory," IEEE Trans. 
Inform. Theory, vol. 44, pp. 2724-2742, October 1998. Commemorative Issue, 
1948-1998, incorporated herein by reference), Rieffel and Polak (E. Rieffel and 
W. Polak, "An introduction to quantum computing for non-physicists " 
http://xxx.lanLgov/abs/quant-ph/98090 1 6. 1998, incorporated herein by reference) 
and Steane (A. Steane, "Quantum computing," Reports on Progress in Physics, 
vol 61, pp. 1 17-173, 1998, incorporated herein by reference). 

As our society becomes more modern and more computerized, there is a need 
to electronically transfer more and more data. To eliminate the limitations of 
conventional data transfer methods, there is a need for a method of data transfer 
according to the quantum information theory. 



SUMMARY OF THE INVENTION 



It is, therefore, an object of the present invention to provide a structure and 
method for encoding/decoding a block of quantum data including removing 
trailing eigenstates from the block that have eigenvalues below a predetermined 
limit to retain leading eigenstates that have eigenvalues above the predetermined 
limit, and encoding the remaining quantum bits retained in the block. The 
remaining quantum bits can also include a linear superposition of the leading 
eigenstates. The predetermined limit is based upon a density matrix of the block. 
This method of encoding produces encoded quantum bits and can further include 
decoding the encoded quantum bits by reversing the encoding. The decoding 
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reproduces the remaining quantum bits and the encoding completely erases the 
remaining quantum bits. Further, the invention can include outputting only an 
encoded or decoded result. 

The invention may further include a method for block compression of 
quantum information which may include projecting a block quantum state into a 
typical subspace having a plurality of eigenstates, encoding the subspace using an 
encoder, and decoding the subspace using a decoder which generates the block 
quantum state using a quantum memoryless Bernoulli source. Projecting of the 
block quantum state into the typical subspace may include analyzing a plurality of 
eigenvalues contained in a density matrix associated with the block quantum state, 
determining a plurality of largest eigenvalues, spanning the subspace wherein the 
eigenstates are associated with the largest eigenvalues to produce a spanned 
subspace, and projecting the block quantum state into the spanned subspace to 
produce a project bloack quantum state that lies in a low dimensional typical 
subspace. The encoder and decoder are quantum-mechanical inverses of each 
other and the decoding is achieved by performing the encoding in reverse. The 
encoding can also include using a fixed-rate quantum Shannon-Fano code to 
compress the projected block quantum state, wherein compression occurs at a per 
symbol code rate that is slightly higher than a von Neumann entropy limit. Also, 
the encoding can include creating a representation quantum Shannon-Fano code 
as a plurality of quantum arithmetic codes and can use the plurality of quantum 
arithmetic codes to compress the subspace containing the projected block 
quantum state. 
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BRIEF DESCRIPTION OF THE DRAWINGS 



The foregoing and other objects, aspects and advantages will be better 
understood from the following detailed description of a preferred embodiment of 
the invention with reference to the drawings, in which: 

Figure 1 is a schematic diagram of reversible circuits; 

Figure 2 is a flow diagram illustrating a preferred algorithm for encoding 
data according to the invention; 

Figure 3 is a quantum circuit implementing the flow diagram of Figure 2; 

Figure 4 is a schematic diagram of a quantum algorithm multiply in which 
i is a classical parameter and B may also be classical, and a schematic circuit 
symbol for the algorithm; 

Figurere 5 is a schematic diagram of a quantum circuit implementing the 
quantun multiplication algorithm of Figure 4; 

Figure 6 is a schematic diagram of algorithms "E" and "D"; 

Figure 7 is a schematic diagram of a quantum circuit implementing the 
block encoder algorithm E in Figure 6 and a schematic symbol for the circuit; 

Figure 8 is a schematic diagram of a quantum circuit implementing the 
block decoder algorithm D in Figure 6 and a schematic symbol for the circuit; 

Figure 9 is a schematic diagram of a quantum circuit implementing the 
Shannon-Fano encoder, the corresponding decoder being obtained by running the 
circuit in reverse; and 
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Figure 10 is a schematic diagram of a hardware embodiment of the present 
invention. 



DETAILED DESCRIPTION OF PREFERRED 
EMBODIMENTS OF THE INVENTION 

The statistics underlying a quantum memoryless Bernoulli source is 
completely captured by the source's density matrix. The fundamental idea behind 
quantum data compression is to analyze the eigen-structure of the joint density 
matrix associated with a block quantum state emitted by the quantum memoryless 
Bernoulli source. 

The following contributions are discussed in greater detail below. The 
invention's first contribution is a quantum-mechanical system for projecting the 
block quantum state onto the subspace spanned by the most important eigenstates 
of the joint density matrix, that is, the eigenstates corresponding to the largest 
eigenvalues. 

In other words, the invention removes the having eigenstates below a 
pre-specified minimum. The remaining eigenstates are the most important and are 
retained. In one embodiment, the invention evaluates the quantum bits by 
determining whether they are typical. The invention computes, in parallel, an 
indicator function that is 0 if the eigenstate is typical and 1 otherwise. 

The determination of what is typical is made based upon the density 
matrix. By making a measurement on the quantum bit associated with the 
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indicator function, with very high probability, the invention projects the block 
quantum state onto the typical subspace spanned by the leading eigenstates. The 
invention represents a strengthening of previous results in that the invention holds 
for fixed block sizes and then delivers a rate of convergence. 

The projection onto the typical subspace wipes out the trailing eigenstates, 
and, hence, the projected quantum state lies in the low-dimensional typical 
subspace. Thus, the invention reduces the number of different dimensions that 
must be compressed, making the processing faster and more efficient. 
Consequently, each leading eigenstate can be represented using roughly the 
algorithm of the dimension of the typical subspace. 

The central problem of quantum data compression is to efficiently 
compute such low-dimensional representations. The invention also includes a 
quantum version of the classical Shannon-Fano code to represent and, hence, 
compress the projected block quantum state using a per symbol code rate that is 
slightly higher than the von Neumann entropy limit. Conceptually, the invention 
achieves data compression by dimensionality reduction. 

As another contribution, the invention proposes arithmetic codes to 
efficiently implement quantum Shannon-Fano codes. The invention's arithmetic 
encoder/decoder uses a certain finite-precision arithmetic process that is inspired 
by classical arithmetic coding. One point of novelty of the inventive quantum 
arithmetic coding is the implementation of a finite-precision arithmetic processes 
in a quantum mechanically reversible fashion. The invention's arithmetic 
encoder/decoder has a cubic circuit and a cubic computational complexity in the 
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block size. The proposed encoder and decoder are quantum-mechanical inverses 
of each other, and constitute a very satisfying example of reversible quantum 
computation. 

The first step of the present invention begins by reviewing the definitions 
of quantum sources and quantum states relevant to the present coding problem. 
Further, the invention presents precise quantum counterparts for the classical 
notions of fidelity and entropy, and describes how encoding and decoding is done 
using quantum computation. 

A classical memoryless Bernoulli source emits a sequence of independent 
and identically distributed symbols each of which is 0 with probability p or 1 with 
probability 1 - p, where 0 < p < 1. The problem of classical noiseless data 
compression is the transmitting of sequences of samples emitted by such a source 
using a minimal number of bits see Shannon, C.E. Shannon, "A mathematical 
theory of communication," Bell System Technical 1, vol 27, pp. 379-423, 1948, 
(incorporated herein by reference) which established that on average each symbol 
can be transmitted in (slightly larger than) H (p) = - p log p - (1 - p) log (1 - p) bits 
with high probability of correct reception, where H (p) is known as the Shannon 
entropy. 

A pure two- dimensional quantum state is known as a quantum bit or 
qubit. The quantum state of a qubit is mathematically represented by a unit norm 
vector in a two - dimensional complex vector space (called a Hilbert space) 
written as H 2 . A qubit may be thought of as a column vector, and is usually 

written using Dirac's ket notation; for example, ) denotes a qubit. The 
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conjugate transpose of ^ , namely, ^ , is written in Dirac's braket notation as 

x . The inner product between an ordered pair of qubits, , is written in 

Dirac's braket notation as ^ |(p ) . The invention writes the fidelity between a pair 

ofqubits, (4» 9 q>) » as <4»|q> ><<P |* > = |<4»|<P >l 2 . Letfoo) and |^> denotetwo 
arbitrary qubits. A quantum memoryless Bernoulli source emits a sequence of 
independent and identically distributed symbols each of which is \ty o ) with 

probability ofp or j<|> 1 ) with probability of 1 - p, where 0 < p < L The per-symbol 
distribution of this source is described by the density matrix: 

P = P\h X A| + 0 * P)\A )(A\ , where denotes the 2 x 2 matrix given 

by the outer product between the vector |(j> > and its conjugate transpose (§\ . 

The problem of (pure state) quantum noiseless data compression is the 
transmitting of such sequences of symbols with high fidelity, using a minimal 
number of quantum bits. According to Schumacher's theorem, (B. Schumacher, 
"Quantum coding," Physical Review A, vol. 51, pp. 2738-2747, 1995, [1 1], 
incorporated herein by reference), on average each symbol can be transmitted in 

(slightly larger than) s(p) = -Tr(p logp) quantum bits with high probability of 

correct reception, where s(p) is known as the von Neumann entropy. A 

surprising contrast between the classical and the quantum cases is that 
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S{p)< H(/>) ^ w here the equality is achieved if and only if the quantum states 

\$ o) and 1^ are orthogonal. Intuitively, this holds since two non-orthogonal 
qubits cannot be distinguished with certainty by measurement. The invention will 
let P and E denote the probability and the expectation, respectively, with respect 
to the quantum memoryless source. 

The invention shall focus on compressing a block of T| symbols emitted 



by the quantum source. Let 



V [M] > s V 1 > ® | V 2 >® • • • Y „ > be a sequence of 



symbols emitted by the quantum memoryless source, where ® denotes the tensor 
product, and \V\ ) represents the rth sample from the source, a random state which 



is either 



^ 0 ^ with probability p or <|> i ) with probability 1 - p. 
This notation is slightly unconventional. Usually, a quantum state written 

using the ket notation \) is definite or pure: it has a von Neumann entropy of 
zero. However, as shown above, the invention shall find it convenient to use 
similar notation to denote mixed or random quantum states (random mixtures of 
pure states); such states will be written with their label in bold, in analogy to the 
often used classical notation for random variables. 

Let x denote a binary string of length 11 . The invention will let ^ denote 
bits /* through j of X. As an alternative to xu,k] , the invention may sometimes 
write xik]to represent the first k significant bits of X. Similar notation, when used 
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with qubits, should be clear by analogy. For example, the invention writes 
^ as ^ [v] ^ . For every a > 0, the invention writes ^ a — ^ and 

a times 

U = 11 - l 

"Vtimes^ . When it is clear how many zeros or ones are necessary, sometimes 
the subscript will be suppressed. 

Although |(|> o ) and |<|> i ) may be non-orthogonal, a basis always exists for 
H 2 in which the same per-symbol distribution source that emits one of two 

orthogonal states |0> and |l> with probabilities ^ ° and ^ 1 (instead of ^ and 

1^ with probabilities p and 1 - /?). The special basis j 0),|l)| is defined as 

follows. The density matrix p characterizes the per-symbol distribution of the 
original quantum Bernoulli source and is self-adjoint, positive-definite, and has 
unit trace. Hence, its eigenvalues, say, ^ 0 > % x , are real and nonnegative, and 

sum to one. If A ° Al , choose the states 0) and 1) to be eigenvectors of p 
corresponding to A, 0 and X \ , respectively. By construction, these states are 
orthonormal, and form a basis for H 2 . If A, 0 = X \ , the select j|o, 1)| to be any 

orthonormal basis of H 2 . In this basis, the original symbols \§ o) and \§ i ) are 
givenasj^) =(0^)10) +<l|^)|l> 

, and the density matrix can be written as 
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p= UloXOl+^lXll , where l 0 = \- l l = p [< oj^ 0 >l 2 + (l- p)|<0^,)| 2 • 

Furthermore, the invention can now write S(p) = H(Xo). 
For 1 < / < n, the mixed quantum state |vj/i) can be written as 

M = <0|i|/i>|0> + (11^)11), 
where (0|i|/i) and are random quantities such that 

E KOMI 2 = Xo and £|<l|\i/i>| 2 = X,. (2) 
The sequence of symbols \y w ) in Eq. (1) is a mixed quantum state in the Hilbert 

space H f 7 = ® i = i m . Using properties of the tensor product, the foregoing can 



be written as 



i 

= ®r =1 x <Xii^>iXi> 

x.=o 

= E(X|^)|X), 



where 

1=1 

Now, 

• (2 (a) J? (b) J? 



(3) 



(4) 



where (a) follows from independence and (b) follows from Eq. (2). The In 
quantum states |X>, X e {0, 1}", can be thought of as the eigenstates of the tensor 



AM9-99-0152 13 



product density matrix p m = S U p , and the numbers A(X;X 0 ,h) as the 
corresponding eigenvalues. Note that the eigenstates |X), X e {0,1 }" constitute an 
orthonormal basis for theHilbert space h® 7 . 

It follows from Eq. (3) that the invention writes the message |v|/ [nJ ) to be 
encoded as a linear superposition of the 2" eigenstates (e.g., |X>, X € {0, l} n ). The 
"randomness" of the message is completely contained in the coefficients <X|v]/ [n] >, 
and the eigenstates are not a function of the particular message to be transmitted. 
Physically, the randomness is embedded entirely in the complex amplitude 
associated with each path or eigenstate. 

The encoding and decoding of classical information is specified by a 
mapping between bit-strings. Similarly, for quantum information, the invention 
specifies a mapping between quantum states. However, additional reversibility 
constraints must be satisfied with quantum states. For example, a reversible 
transformation must conserve energy. Since quantum states are mathematically 
represented by vectors with unit norm, reversible transformations must preserve 
the norm. It also turns out that with the appropriate description of the system, the 
most general transformation preserves orthogonality between states. If the 
quantum states in H<> as 2 n - dimensional column vectors, then most general 
transformations are described by 2" x 2" unitary matrices acting on the Hilbert 
space of the quantum states. Again, a unitary matrix is one whose conjugate 
transpose is its inverse. 
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This model of computation subsumes classical computation, because 
mappings between bit-strings can be described as permutation matrices acting on 
the basis elements of the Hilbert space. Unitary transforms are always invertible 
or reversible. All irreversible (classical) computation can be made reversible with 
only a polynomial amount of overhead (see C. H. Bennett, "Logical reversibility 
of computation " IBM J. Res. Dev., vol. 17, pp. 525-532, 1973, incorporated 
herein by reference). However, not all unitary transforms represent reversible 
classical computation. Not all unitary transforms can be described by permutation 
matrices. A unitary transform can be completely specified by its action on all the 
basis elements of a Hilbert space. Transformations which are not permutations 
take basis elements to superpositions of basis elements; these are at the heart of 
the speedup of quantum computation and quantum error correction. 

Quantum algorithms are generally very difficult to construct, but choosing 
the eigenstates |X>, X € {0,1 as the basis vastly simplifies the descriptions of 
the inventions encoding and decoding transforms. In this special basis, the 
invention only employs unitary transformations which are permutations of basis 
elements to achieve the inventions goal. These transforms shall be applied to 
input states which are generally in non-classical superpositions of basis elements. 
As suggested by Deutsch, (D. Deutsch, "Quantum theory, the Church-Turing 
principle and the universal quantum computer," Proc. R. Soc. Land. A, col. 400, 
pp. 97-1 17, 1985, incorporated herein by reference), it is convenient to think of 
what happens as being "quantum parallelism". For an input |<j>> = a|0)+b|l), a 
computation U produces U \§) = aC710>+bC^l), by linearity. For example, two 
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"classical" computation happen in parallel, one with input |0> and the other with, 
|1>, with the two computational paths being weighted with complex amplitudes a 
and b, respectively. Similar observations hold for arbitrarily large states. As long 
as U is simply a permutation, these different paths never interfere and a coherent 
quantum state is maintained. 

The invention symbolically describes encoding and decoding unitary 
transforms for quantum information using algorithms which at first glance look 
very classical, but in reality, are specially constructed to be quantum. Three 
characteristics make the invention's algorithms quantum-mechanical. First, they 
are reversible, which is required as previously explained. Second, they completely 
erase their inputs, which is a necessity because quantum states cannot be cloned 
(see W. K Wootters and W.H. Zurek, "A single quantum cannot be cloned," 
Nature, vol. 299, pp. 802-3, 1982, incorporated herein by reference) and D. Dieks, 
"Communication by EPR devices," Physics Letters A vol. 92 (6), pp. 271-272, 
1982, incorporated herein by reference), and thus there is no sense to a sender 
sending a faithfully encoded quantum state elsewhere without erasing their own 
knowledge of that state in the process. Third, the invention produces no 
information other than the encoded (or decoded) state, which allows 
differentiation between computational paths. Producing such entanglement would 
ruin the superposition which is being encoded, because any potential for obtaining 
"which path" information implies the existence of a physical measurement which* 
would (at least partially) collapse the superposition state. Fundamentally, this non- 
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disturbance requirement is deeply related to the no-cloning theorem, and it is a 
subtle, but very important point. 

The invention employs for clarity of exposition, quantum circuits, which 
succinctly capture the same information as the algorithms, and often effectively 
convey additional structural information about the procedure. A wide body of 
knowledge about quantum circuits exists (see, A. Barenco, C.H. Bennett, R. 
Cleve, D.P. DiVicenzo, N. Margolus, P.W. Shor, T. Sleator, J. Smolin, and H. 
Weinfurther, "Elementary gates for quantum computation," Physical Review A, 
vol. 52, no. 5, pp. 3457-3467, 1995, incorporated herein by reference, and A. 
Barenco, "A universal two-bit gate for quantum computation," Proc. R. Soc. 
London A. Vol. 449, no. 1937, pp. 679-683, 1995, incorporated herein by 
reference), but only the subset shall be drawn from it, which is convenient for 
describing reversible classical circuits, including the controlled-NOT and swap 
gates, as shown in Figure 1. 

Figure 1 illustrates two quantum gates which are used in later quantum 
circuits. The first gate is the controlled-NOT gate, which produces x' = x and y' = 
© y (© denoting (bitwise) addition module 2). The second gate is the swap gate, 
for which x f = y and y' = x. Time goes from left to right. A final useful notation 
for expressing the coding procedure is that, given a fractional number (e.g., ^,0 < 
£ < 1), let £ = 0.£i 0 • * * , denote a binary representation of the number. This 
number can be associated to a pure quantum state \Q = |£i> ® |£ 2 ) ® ® • • ■ 

in the infinite-dimensional Hilbert space H i . This allows for representation of 

fractional real number as a quantum state. 
AM9-99-0152 17 



As shown in Figure 2, a symbolic system or "psuedocode" for computing 
Eq. (6), discussed below. w <~ " denotes an assignment operation. When 
describing a pre-existing state or comparison operation, the invention uses ' A 
temporary quantum register |w) of length [log n] is used. This register is initialized 
and finalized to |0ri O gn]>. The precise value of x should satisfy Eq. (8), below and 
will be specified further below in Eq. (16). 

Figure 3 illustrates a quantum circuit implementing the system in Figure 2. 
The gates labeled as U + and tl implement lines 2-4 and lines 10-12 of the system, 
respectively. These gates are quantum-mechanical inverses of each other. The gate 
U> implements lines 6-8 of the system. As in "Barenco", supra, the invention 
generally uses rounded symbols to denote the control qubits, and boxed symbols 
to indicate the targets, with the exception of"®" which always sits on a target. 
The / losn1 notation indicates a wire bundle with Tlogn] qubits. This is proven in 
the adaption of the technique in Hoeffding (W. Hoeffding, "Probability 
inequalities for sums of bounded random variables," American Statistical 
Association Journal, vol. 58, pp. 13-30, 1963, incorporated herein by reference). 

The basic idea of quantum data compression is that the eigenstates 
associated with smaller eigenvalues can be discarded without incurring significant 
loss of average fidelity. This goal will be attained by employing a measurement of 
a certain quantum observable associated with the given message, described further 
below. 
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Let w(X), X € {0,1 }\ denote the Hamming weight of the string X, for 
example, the number of ones in the string. It follows from Eq. (4) that it can be 
written 

A(X;^,)= Al HX) A? iX) (5) 
Since k 0 > X u it follows from Eq. (5) that smaller the Hamming weight of an 
eigenstates the larger the eigenvalue associated with the eigenstate. Let x > 0 
denote a truncation threshold. Let G T and B T denote the sets of "good" and "bad" 
eigenstates such that 

G r = {X|w(X) < r} 
B r = {X|w(X) > r}. 

With appropriate values of x, the subspace spanned by the good eigenstates, 
namely, 

span{|X>|X g G T } 

becomes the typical subspace that contains most of the information present in an 
average quantum message. 

For every eigenstate |X>, X e {0,1 } n , let I T (X) denote the good-bad 
indicator function such that 

fOifX € Gx, 
Il(X) " (l ifX « B,. 

The following transformation is now computed: 

|X,0> - |X,I t (X)>. (6) 
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The invention exhibits a quantum system for computing Eq. (6) in Figure 
2, which is implemented by the quantum circuit in Figure 3. This system makes 
use of subroutines previously described in the literature (see R. Cleve and D.P. 
DiVicenzo, "Schumacher's quantum data compression as a quantum 
computation," Physical Review A, vol. 54, pp. 2636-2650, October 1996, 
incorporated herein by reference) for conditional addition and subtraction, and 
comparison. Using Eq. (3), the action of the system on the quantum message can 
be written as 

kw.o) -> E<x|r[*]>|x,ir(X)> ■ \fm,h), 

Xe{0,l}i ^ 

where 5 is an output state in which I T is now a function of ^ ^ and thus, 
in general, is entangled with it. Let 

m m 

{|Ir> =|0» or =|1)} 

denote the two possible events corresponding to measuring |I T ) to be |0> or |1>, 
respectively. The truncation threshold t is determined to ensure that the 

m 

probability of the event ~ ^ is close to 1. 

Assuming that y 2 <Xo< l,fora fixed n > 1 and a fixed 8 > 0, it is set that 



T > 



8 ^ 



V \ogAo/A\J 



(8) 



then 



= |0» = l-P{|Ir> = |1» > l-Z^o'^ . (9) 
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The invention adapts the technique in Hoeffding (29, W. Hoeffding, "Probability 
inequalities for sums of bounded random variables," American Statistical 
Association Journal, vol. 58, pp.13-30, 1963.) incorporated herein by reference: 



r{|ir|i» = I„ B ,E<*|r w ) 



(b) 



(c) 
< 



mm 



y > 0 



£ 2^a(xaaM^M)) A (,;4 ) ,4 i ) 

(x|-k)gA(XUo.-*i)i9(s(/>)+*)} 
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<min 



y >0 



= mm l. 



7 

2-^f| (/l 0 2''(- lo 8^- s W) + /l^^- 10 ^--^))) 



Y>0 



(d) 



<min 



^ > 0 



lt(2< 1/8 )r 2 (iog(^a,)) 2 ) 



(e) 2- 2 " y2 '('°8(- i o'' ! ')) 2 



where (a) follows from Eq. (4); (b) follows since o ,0<9< q is a 
decreasing function of 0, and, hence, by using Eq. (8) 



X\ XJ~ T < ^i+^OogUo/zii))^ i-nM-nsiiogUo/h)) 
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1 



— X \ X\ 



o 



v v 



7 <y/(iog(v*)) 



_ 2""(s(pH 



(c) holds for all Y> 0, since 



2>'(-logA( 7 ,.! 0 J 1 )- 7 (s(,>M) 



(d) follows from [W. Hoeffding, "Probability inequalities for sums of bounded 
random variables," American Statistical Association Journal, vol. 58, pp. 13-30, 
1963, incorporated herein by reference, (4.16)], if 1 / 2 < A 0 < 1 ; and (e) follows 
by selecting the minimizing value 



AS 

r = 



(logfe/i,)) 



2. 



10 Observe that \f [„]) and | I r ) in Eq. (7) are, in general, entangled. Hence, 

a measurement on the last qubit will irreversibly effect the first n qubits. Precisely, 
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using yon Neumann's postulate ( J. Von Neumann, Mathematical Foundations of 
Quantum Mechanics. Princeton, USA: Princeton University Press, 1955. 

Chapter VI), incorporated herein by reference, the effect on 1^ [ >?] ) of measuring 



is the following: 



*«'{z\r[,))\ 2 



-Jz **** 



(z 



In other words, if the event 



i,} a |o> 



occurs, then W[n)j will collapse to the 



renormalized projection of the message onto the subspace spanned by the 



good eigenstates, otherwise f\ n \] will collapse to the renormalized projection of 

the message ¥[ n ]j onto the subspace spanned by the bad eigenstates. Thus, with 
high probability, the bad eigenstates are discarded. 

It follows from Theorem 3.1 that the event j| I r ) S 1 0) j occurs with very high 
probability. When this event occurs, the invention now illustrates that, the 

collapsed state ^]) is not much different from the original message 
that is, the average fidelity between the two is close to the maximum possible 
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value of 1. Recall that the average fidelity is the probability that the message 

$ [ n \ ) passes a test for being the same as the original message | ^[7]) , whence 

the test is conducted by someone who knows the original message, see, 
Schumacher (1 1, B. Schumacher, "Quantum coding," Physical Review A, vol. 51, 
pp. 2738-2747, 1995.) incorporated herein by reference. 

If the previous hypotheses holds, then 



^wkw>i 2 i{|i.> 2 io>} 



>\-2 (togv*) 2 



This is shown in the following. 



<^wK]>l 2 |0r>a|0» 



= E 



re{o,i}*<Mo,ir X J/ 1 K 



(a)E 



*«{o,i}' 
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= E I (z\%])\ 2 



(d) i- 2>(*;4a) 



> 1_ 2 (log/ioAii) 2 



where (a) follows by using the orthonormality of the eigenstates; (b) follows from 
Eq. (10); (c) follows from Eq. (4); (d) follows by applying the binomial theorem 
t0 1 = (i 0 + ii) 7 ; md ( e ) fol l° ws fr° m tne foregoing. 

The foregoing, represents a strengthening of Schumacher's pioneering 
result in that they hold for fixed block sizes and they deliver a rate of 
convergence. 

The invention proposes the following scheme for transmitting the quantum 
message 
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compute Eq. (7) 



measure 



» l p<H°>}) 



then 



transmit 
else 
do nothing 

It follows that the above scheme has high average fidelity with high probability, 
and only an exponentially small probability of failing to transmit any information. 



This can be explained as follows: the desirable event 



lr)*|0>}) 



occurs with 



probability close to 1. And, in the case of this even, only the bad eigenstates are 
discarded. 



From now on, it is assumed that the event 



I r )s|0)} has occurred, and the 



following focuses on transmitting f^). It follows from Eq. (10) that ^ ^ lies 

in the typical subspace spanned by the good eigenstates. The invention selects the 
truncation threshold such that the typical subspace has dimension at most 
2>i{s{p)+s)+i which is much less than the original dimension of 2rj. Hence, by 
appropriately "relabeling" the leading eigenstates, the invention should be able to 
represent, and, hence, compress the n qubit message to ?j(s(p) + s) + 1 
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qubits. The main problem, which is now tackled, is how to compute such a 
dimensionality reducing or relabeling transformation efficiently. 

The eigenvalues Xq and X\ are real numbers, and, when represented as 
fractional binary numbers, may require an infinite precision to represent. Since, in 
practice, the invention can only store and manipulate a finite number of bits, from 
now on, the invention approximates the eigenvalues using fractional numbers with 

q significant bits after the binary point. In particular, the invention lets % \ 

denote the fractional number obtained by truncating all but the q most significant 

o / <f 

bits of X 0 . And, the invention lets Xi = X\ + \Xq-Xo 



o 

Since, X Q + X\ - 1 ? it follows that A l has at most q nonzero significant bits, and 

the remaining bits must be zeroes. Furthermore, 

o o 
X o+ X i = 1. 

In the rest of the disclosure, instead of the original eigenvalues, the invention will 

use X o m & X x • To be sure, such an approximation will slightly increase the per 
symbol rate needed for compression by 



D 



•o 

V 



X o I = A 0 logl X 0 1 X o J + X l logl X { I X i 
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The quantity D(*||«)is known as the relative entropy or as the Kullback-Leibler 
distance. This increase in the per-symbol rate can be made as small as desired by 
selecting a large enough q. However, the invention will subsequently demonstrate 
that the amount of quantum hardware required to implement the encoders and 
decoders will increase quadratically in q. 

The invention now introduces a quantum "encoder* transformation that 

transforms each eigenstate z)>Z e {°> as follows: 



x, o^|o,,c( z )) 



(12) 



where, for a > 0, 0 a represents a string of a zeroes, and 



HwY^z (13) 



( 0 oN l 

where A \%\A o,/t 1 is obtained from Eq. (4) and -< denotes some total order 



on the strings in {0, 1 }". The invention will specify a computationally 
simple-to-implement lexicographical order below. Observe that for every 

eigenstate | z), c(z) is a number in the real interval (0, 1 ). Hence, given C<», the 

t 

invention writes |c(/)} using the terminology of Section 2.5. Intuitively, C( x ) is 

the sum of the eigenvalues of all eigenstates of length n that are less than or equal 
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to the x in the total order -< . Since C(X) is a monotonically increasing function of 
the eigenstates arranged in lexicographical order, it is uniquely decodable. In other 
words, the transformation 



(14) 



exists for every eigenstate z ), z e {o, l}". Hence, Eq. (12) is reversible, and 

can be implemented as an unitary transformation. 

Each eigenvalue in the sum Eq. (13) is a product of n numbers each of 
which has a precision of q bits. Hence, each eigenvalue can be written as a 
fractional binary number with at most nq nonzero significant bits. Finally, this 
implies that, for each eigenstate, the number C( x ) has precision no more than nq 

bits. In other words, the encoder is a unitary transformation from H 2 to a 

2"-dimensional subspace of H 2 . However, since generally q > 1, this hardly 
constitutes data compression. The invention now achieves compression by 
truncating a large number of nonsignificant bits of C(x). 

For a given truncation parameter k > 0 and a given eigenstate 

LA y g / Q j\' 

' i ' i the invention defines the truncated encoder transform as 
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Eq. (15) 



where C{x\k] denotes the truncation of C(x) to the k most significant qubits. 
Observe that only the first k qubits on the right-hand side depend upon the 
eigenstates, and, hence, only these bits need be transmitted. All the information 
present in the eigenstate x has been captured into C(x) in a reversible fashion. 
Hence, compression does not cause any loss of information. 

Consequently, the encoder in Eq. (15) maps messages of a fixed-length n 
to codewords of fixed-length h In other words, the encoder is a unitary 

transformation fromH 2 to a subspace of H 2 • 

The decodability of the untruncated map in Eq. (12) is immediate from the 
fact that Ox) is a monotonically increasing function of the eigenstates arranged in 
the lexicographical order. In contrast, the decodability of the truncated map in Eq. 
(15) is a delicate matter. If k < n, then the truncated map cannot hope to correctly 
decode all the eigenstates. However, fortunately, the invention only needs to 
correctly decode the good eigenstates. The invention can discard the bad 



eigenstates, since we have assumed that the event 



lr> = |0>}) 



has occurred, and 



in this case, bad eigenstates have already been discarded. So, there was no need to' 
either encode or decode the bad eigenstates. 
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The invention now establishes that if the threshold parameter T and 
truncation parameter k are carefully selected, then inverse of Eq. (15) exists for all 
the good eigenstates. 
If, 



T = 



7 



X\ + 



logio/zi,. 



+ 7<?+log(V/0- 



(16) 
(17) 



Then, there exists a decoder such that, for every X in G T , 



(18) 



This is shown in the following. Given the encoding C( l w _ k ) 0 f the 



eigenstate \x) , the invention defines the corresponding decoded or reconstructed 



eigenstate as 



j that satisfies the following two inequalities: 



Eq. (19) 



i=k+l 



C W[k] + t 2- i <c(f) + A(f;/ 0 ,i') 

i=k+l v J 



Eq. (20) , 
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In general, owing to truncation, the decoded eigenstate % need not equal the 
original eigenstate X. The invention now shows that for values of f as in Eq. 
(16), for values of A: as in Eq. (17), and for all good eigenstates, the inequalities 
Eq. (19) and Eq. (20) are satisfied if and only if % = X ■ Suppose that X = X 
In this case, the first inequality Eq. (19) is trivial, and holds for allXin {0,1}". 
Now, observe that the second inequality Eq. (20) holds if 

A[7;io,iiJ = A^;io,/lil>2-k> £ 2-«. 

i=k+l 



It follows from Eq. (4) that the second inequality Eq. (20) holds if 

( o\ Ax)/ o\ iMx) 

Ui Uo >2- k . 



The invention would like the above inequality to hold for all good eigenstates. 

[a i) [to) ,o<e< tj, 

Since , is a decreasing function of 0 , it is sufficient 

that the above inequality holds for the good eigenstates corresponding to the 
smallest good eigenvalue. If the invention selects T as in Eq. (16), then the 

smallest eigenvalue is larger than y X lj y A o | 
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Hence, the invention requires that I ^ 1 J I ^ °l ^ 2 _k . 



Equivalently, the invention requires that 



WW 



f *\ 



( o\ 
X o 



Tj-r 



>2-k. 



It follows from simple algebraic manipulations that the above holds if 



10 



k> 7S(/?)+7D 4 



/U U^+log^o/zi,). 



This is exactly the requirement in Eq. (17). 

The invention now establishes the converse, that is, if X * X , then both 
Eq. (19) and Eq. (20) do not hold. There are two cases: either / -< % or X ~*> X • 
In the former case, Eq. (19) cannot hold and in the latter case Eq. (20) cannot 



hold. 



Observe that the desired encoder transform in Eq. (15) annihilates the 



quantum state X. This is necessary since both | X ) and ^W[ k ] ) contain the 
same information, and since quantum states cannot be cloned, it is impossible to 
faithfully transmit weighted superpositions of different \x) without the sender 
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obliterating her knowledge about it in the process of transforming the state into a 



weighted superposition of 



c(z) 



[k] 



Observe that the untruncated map in Eq. (12) and the truncated map in Eq. (15) 
map one eigenstate to one encoded state. Hence, they can be thought of as unitary 
transforms that are permutations of the basis states. 

So far, the invention has specified the desired encoder Eq. (15) and the 
corresponding decoder in Eq. (18) in terms of the eigenstates alone. For the sake 
of completeness, by using linearity of the encoder and the decoder, the invention 
now describes their action on the quantum message of interest: 



encode 

-> 



I (zi b ))o n ,c(z) fkl i 



[k] 1 W-k> 



transmit 



prepare 



I <zr[,))o,,c(z) [k] i n - k ) 



ZeGr 



decode __, / i 
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where the invention has implicitly used the fact that a measurement on the qubit 



of equations capture the essence of the invention. The eigenstates are encoded, 
transmitted, and decoded without any loss of information. Conceptually, encoding 
and decoding are reversible operations that are conceptual inverses of each other, 
hence, no information is lost in the entire process. 

The invention now proposes quantum algorithms and associated quantum 
circuits to efficiently realize the encoder in Eq. (1 5) and the corresponding 
decoder in Eq. (18). 

First, the invention considers the computation of the function C(x) in Eq. 
(13). A straightforward algorithm for computing C(x) by explicitly performing the 
summation would require an exponential amount of complexity in the block size 
n. One of the main contributions of classical arithmetic coding is to observe that if 
the invention selects the total order -< in Eq. (13) to be the following 
lexicographical order, then the function C(x) can be efficiently computed. If 

t= $ ,and Z s Xi%2-..Zn are in {0,1}" then the invention says that 



Under this definition of the total order "< the invention can write the function 
C(x) recursively as follows, see, [J. Rissanen et al., supra, (1)]: 



I r ) has been made, and the event 



{Mlo}} 



has occurred. The foregoing series 



ifandonlyif L £ 2W < S X^~^ 



i=l i=l 



20 
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for i=l to n do 

if(JB = 0) 

o 

C(X) = C(X) x Ao 
else 

endif 
endfor 

Instead of the lexicographical order in Eq. (21), the invention can also use 
the following dual order. If £ s - . ^ and % = %\%r ~ % n are in {0, 1 } n , then 
the invention says that 

dual 7} rj 

£ •< jifandonlyif £ £ 2 -i + i < £ ^2" i+1 . 

i=l i=l 

Under this dual definition, the invention can also write the function C(x) 
recursively, see, (Rissanen et al, supra, (2)). Although both the recursions are 
amenable to a quantum implementation, the recursion corresponding to the total 
order in Eq. (21) turns out to slightly simpler and, hence, is used in this disclosure. 

Important parts of the encoding and decoding algorithms are 
multiplication and division, respectively, and in order to build the quantum 
coders, the invention must first construct quantum algorithms for such arithmetic. 
Suitable addition and subtraction circuits have already been described in the 
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literature (D. Beckman, A.N. Chari, S. Devabhaktuni, and J. Preskill, "Efficient 
networks for quantum factoring," Physical Review A, vol. 54, no. 2, pp. 
1034-1063, 1996, http://xxx.lanl.gov/abs/quant-ph/9602016, and R. Cleve and 
D.P. DiVincenzo, "Schumacher's quantum data compression as a quantum 
computation," Physical Review A, vol. 54, pp. 2636-2650, October 1996, each 
incorporated herein by reference), but appropriate multiplication and division 
algorithms have not been. These are described below. 

The invention presents in Figure 4 an algorithm to multiply 

| A ) ? B ), R), i) that takes the following inputs: (a) a fixed index i, / = 1, 2, • • • 

,n 9 (b) nq qubit register A } such that all but the first (M) q qubits are zeroes, (c) 

q qubit register |B ) 9 and (d) q qubit register | R) . The algorithm also requires a 

nq qubit temporary register | y) that is initialized and finalized to 0^ ). The 



algorithm computes 



A,R>-> AB + 2-0-»q-iRn >, 



where multiplications and additions are to be interpreted by treating A, B, and R as 
fractional binary numbers. A quantum circuit which implements the algorithm is 
shown in Figure 5. The steps in Figure 4 correspond simply to a quantum version 
of the classical algorithm for multiplying two numbers A and B where there is an 
additional quantity that is added at the end. 

« 

The invention terms the conjugate inverse of this algorithm as divide 
a),|b),|r) , i . Given a nq qubit register U \ such that all but the first iq 
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qubits are zeroes, a q qubit register g ) , and a q qubit register |r) that is 



initialized to 



0 q ) , the circuit divide ( A ) , |b ) , |R) , ij uses a nq qubit temporary 



register | T) that is initialized and finalized to |o 7q > and divides A by B up to the 

first (/ - 1 )q bits, and stores the quotient also in A, and keeps the q qubit remainder 
mR. 

The invention now uses the ideas from arithmetic recursions, and the 
above circuits for multiplication and division to construct building blocks for the 
desired encoder in Eq. (15). In Figure 6, the invention presents two recursive 
algorithms "E" and "D." Formally and literally, these algorithms are inverses of 
each other: lines E2-E8 are literal inverses of lines D7-D13, lines E9-E13 are 
literal inverses of lines D2-D6, and, finally, the for loop in the algorithm E 
processes the message symbols in the original order from 1 to n while the for loop 
in the algorithm D emits the message symbols in the inverse order from n to 1. 
The invention exhibits quantum circuits for implementing the algorithms E and D 
in Figures 7 and 8, respectively. Observe that these circuits are also 
quantum-mechanical inverses of each other. 



The invention intends to use the algorithms E and D with two different 
sets of inputs. The invention now explains the functionality of these algorithms on 

« 

the first set of inputs. 
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Let \z)> X e ' denote m y ei g enstate - The algorithms D and E, 



respectively, compute the following maps: 



D,:0„c(ir),O w >-> U,o„,o„>, 



Eq. (23) 



Eq. (24) 



With the inputs as above, Ei is a quantum version of the arithmetic recursion 
presented above. The desired assertion for Di follows by observing that it is a 



This furnishes a way of implementing Eq. (12) and its inverse Eq. (14). 
Recall, however, that to achieve compression the invention is interested in 
implementing Eq. (15). The obvious strategy of first implementing Eq. (12) and 



simply transmitting the k most significant qubits of C( %) ) does not work, since 



these k qubits are entangled with the nq-k least significant qubits of C( % 



Hence, a measurement on these nq-k least significant qubits will irreversibly 
change the k most significant qubits. To avoid such an accident, the invention 
must erase the nq-k qubits. This is the central difficulty that the invention must ' 
overcome. The invention now explains the functionality of the algorithms E and D 
on the second set of inputs. 
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literal inverse of Ei . In both of these cases, the quantum register | R X R 2 * • • R n 



always remains in the same initial state 0 






event 



Suppose that a measurement on the qubit 1 1 r ) has been made, and the 
1 1 r )~ 1 0)} has occurred. Let I ^ )' ^ 6 ® r , denote any good eigenstate. 



The algorithms D and E, respectively, compute the following maps: 



A:0, f ^) W V*,° W ) >\z,K,R&-R n ) (25) 



£ 2 :|j, 0 W , R X R 2 -R n ) 



(26) 



The invention establishes the assertion for D 2 in detail. The desired 
assertion for E 2 follows by observing that it is a literal inverse of D 2 . Fix a good 



correctly; the gist of what follows is that not only may | % ) be decoded correctly, 
in fact, it may be decoded correctly in a sequential or recursive fashion. For a 
index i, i=l, 2,..-, ^recall that 1^)= \ ZlZ2 .- Zi ). 

The invention first shows that lines D7-D13 in Figure 6 behave as desired. For 
any / =1 , 2, • • •, r\ and for real number C, observe that if 



eigenstate | %) = \ X\ X2 • • • Zn ) • As shown above, I % 



can be decoded 





(27) 
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then, after the computation in lines D7-D13, the invention has 

c(j [M] )< C< c(j [M] )+ A( / r M1 ;io ? i i) 



However, the invention has from Eq. (19) and Eq. (20) that, if the invention sets 

C = C( x) [k] + X Tk+i ' ^ en * e ine( l ualit y Ec *' ( 27 ) holds for f " n ' Hence > 
by induction, the inequality Eq. (27) holds for all i=\ , 2, • * % ?/. 
Step 2: The invention now shows that lines D2-D6 behave as desired. For any i-l, 
2, • • rf and for real number C, if the inequality Eq. (27) holds, then 

C > X o if and only if c(/ [i] ) > X o (28) 
The "if part of Eq. (28) follows trivially from Eq. (27). To see the "only if 1 part, 

observe that if Q> x\ , then, by Eq. (27),c(/ [i] )+ A^ji o,X\)> Xo . 



C (l[i]) > ^°- A (l[i]^ 0 ^ 1 

Hence, 



Hence, again by Eq. (27), 



The only allowed values of c( J[*j ) ^ at sat * s fy above inequality are 

C(z[i])^ X o as desired. 
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Observe that Eq. (25) is almost the desired decoder Eq. (18) except for the 

"remainder" |RiR 2 -R 7 ) , which is left over. Once again, the decoded state) x) * s 

entangled with this remainder, and, hence, the remainder must be erased. 
Similarly, Eq. (26) is almost the desired encoder Eq. (15) except that it 
requires the above left over remainder as an input. 

It follows from the above discussion that the algorithms described above 
do not, in themselves, yield either the desired encoder Eq. (15) or the decoder Eq. 
(18). The invention now presents an algorithm, in Figure 9, the desired encoder. 
The desired decoder is obtained by literally running the encoder in reverse. 
The circuit in Figure 9 is started by applying the transformation Ej 

E,:|;r,o w ,o„)-> \o n ,c(%),o m ). 

After the k most significant qubits of ^ are copied (of course, they are not 

truly copied in the classical sense, since qubits cannot be cloned; they are 
entangled with an auxiliary set of qubits prepared in the | q) state), the output of 

Ei is acted upon by the transformation Di 

D Y ]o^c( z )fi n }^ |jr,o w ,o TO ). 

This has the effect of annihilating all the nq qubits of j C(/)^ . However, it 
recreates the input quantum state | % ) which must also be erased. Now, by 



employing the k copied qubits 



C(/)j k j ^ , the invention can apply D2 
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D 



The quantum state | produced at the output of D 2 is used to erase the 



same 



quantum state produced at the output of Di . Now, by applying E 2 to the output 
produced by D 2 , the invention has the desired output: 

E 2 :U,0 w ,R 1 R 2 ...R,)-> \%A%\^m-^. 



In the end, the invention is guaranteed that no quantum register in Figure 9 is 

entangled with the final output ^{x ) [k] ) . Hence, the output can now be freely 

transmitted. Observe that the cascade of Ej and Di is the identity map, and, 
similarly, the cascade of D 2 and E 2 is also the identity map. 

The invention now analyzes the complexity of implementing the Ei block 
in Figure 9. The Ei block can be implemented using the circuit presented in Figure 
7. The ">" operator compares a nq qubit register C to a q bit constant. Using the 
TEST-GREATER-THAN circuits [R. Cleve et al., supra], such comparisons can 
be implemented quantum-mechanically in 0(nq) elementary quantum gates. The 
invention has used a "swap" or o operator in circuits for multiply and divide. A 
quantum-mechanical operator that swaps two quantum registers of length q can be 
implemented using 0(q) quantum Fredkin gates (E. Fredkin and T. Toffoli, 
"Conservative logic," International Journal of Theoretical Physics, vol. 21, i 



no. 
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3/4, pp. 219-253, 1982, and H.F. Chau and F. Wilczek, "Simple realization of the 
Fredkin gate using a series of two-body operators " Physical Review Letters, vol. 
75, no. 4, pp. 748-750, 1995, each incorporated herein by reference). For the 
index z, 1 < / < n, the overall circuit for Mi can be implemented in 0(i 2 q 2 ) 
elementary quantum gates. In conclusion, the overall circuit for the Ei block can 
be implemented using 0(n 3 q 2 ) elementary quantum gates. The blocks D u E 2 , and 
D 2 have the same complexity as the block Ei. Hence, the overall encoder in Figure 
9 can also be implemented using 0(n 3 q 2 ) elementary quantum gates. Also, using 
similar reasoning, it follows that the overall encoder in Figure 9 has a 0(r?q 2 ) 
computational complexity. 

While the overall methodology of the invention is described above, the 
invention can be embodied in any number of different types of systems and 
executed in any number of different ways, as would be known by one ordinarily 
skilled in the art. For example, as illustrated in Figure 10, a typical hardware 
configuration of an information handling/computer system in accordance with the 
invention preferably has at least one processor or central processing unit (CPU) 
1000. For example, the central processing unit 1000 could include various 
image/texture processing units, mapping units, weighting units, classification 
units, clustering units, filters, adders, subtractors, comparators, etc. Alternatively, 
as would be known by one ordinarily skilled in the art given this disclosure, 
multiple specialized CPU's (or other similar individual functional units) could 
perform the same processing, mapping, weighting, classifying, clustering, 
filtering, adding, subtracting, comparing, etc. 
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The CPU 1000 is interconnected via a system bus 1001 to a random access 
memory (RAM) 1002, read-only memory (ROM) 1003, input/output (I/O) adapter 
1004 (for connecting peripheral devices such as disk units 1005 and tape drives 
1006 to the bus 1001), communication adapter 1007 (for connecting an 
information handling system to a data processing network) user interface adapter 
1008 (for connecting peripherals 1009-1010 such as a keyboard, mouse, imager, 
microphone, speaker and/or other interface device to the bus 1001), a printer 
101 1, and display adapter 1012 (for connecting the bus 1001 to a display device 
1013). The invention could be implemented using the structure shown in Figure 
10 by including the inventive method, described above, within a computer 
program stored on the storage device 1005. 

The invention has constructed a quantum algorithm for block compression 
of quantum information which is an analog of classical arithmetic coding. In 
contrast to the classical case, the quantum algorithm must take extra care to leave 
behind no residual traces of its past history. The algorithm thus begins by 
projecting the state into the typical subspace, then a sequence of encoding and 
decoding using finite precision arithmetic is done in a manner so as to obliterate 
all possible imprecisions. 

Unlike the classical algorithms for arithmetic coding, the multiplication 
steps used in the invention's algorithm require a linearly increasing precision in 
the block size n. In the classical case, it is known how to implement these 
multiplications using precision that is independent of n [J. Rissanen et al., supra]. 
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It is straightforward to perform this algorithm in parallel, so as to reduce 
the number of time-steps necessary for its circuit implementation. Multiplication 
and addition are known to be in NC(l), and believed to also be in the quantum 
counterpart to this class, so that it is possible to obtain an 0(h) running time 
implementation of the inventions algorithm. Quantum circuits such as the one the 
invention presented may also find use as reversible classical circuits, which 
potentially require much less power for their execution when using technologies 
such as reversible CMOS or charge recovery logic see, (S. Younis and T. Knight, 
"Non dissipative rail drivers for adiabatic circuits," in Proceedings of the 
Sixteenth Conference on Advanced Research in VLSI 1995 (Los Alamitos, CA), 
pp. 404-414, IEEE Comput Soc. Press, 1995), incorporated herein by reference). 
Finally, the invention has considered block arithmetic codes. 

While the invention has been described in terms of preferred 
embodiments, those skilled in the art will recognize that the invention can be 
practiced with modification within the spirit and scope of the appended claims. 
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